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We consider the Ferromagnetic-String-Formation-Probability correlation func- 
tion (FSFP) for the spin-| Heisenberg XXZ chain. We construct a completely 
integrable system of integro-difference equations (IDE), which has the FSFP 
as a t- function. We derive the associated Riemann-Hilbert problem and ob- 
tain the large distance asymptotics of the FSFP correlator in some limiting 
cases. 
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1 Introduction 



In this paper we continue our investigation of a particular zero-temperature correlation 
function of the XXZ Heisenberg model in the critical regime — 1 < A < 1 in an external 
magnetic field. The XXZ hamiltonian is given by 

n £ rrjnj. ; + a]a] . , + A (a*a* +1 - 1) - ha] , (1.1) 

3 

where the sum is over all integers j, o a are Pauli matrices and h is an external magnetic 
field. For later convenience we define A = cos(2r/) = w , where | < r\ < n. The FSFP 
correlation function is defined as follows 

m 

P(m) = {0\ HPj\0), (1.2) 

3=1 

where |0) is the antiferromagnetic ground state and Pj = + 1) is the projection 
operator onto the state with spin up at site number j. The physical meaning of P(m) is 
the probability of finding a ferromagnetic string (i.e. m (adjacent) parallel spins up) in 
the ground state |0) of the model ( |1 . 1| ) for a given value of the magnetic field h. From 
a mathematical point of view this correlator turns out the simplest one to be considered 
(see i). 

In a recent paper || we have derived a representation of the correlation function ( |1.2|) 
as a determinant of a Fredholm integral operator 



P(m] 



(0 |det (l + 7< 



m) 







1.3) 



det (l - IK) 

Here K and t/( m ) are integral operators acting on a function /(A) as 

(yM/) (A) = vW(\\n)f(n) , (Kf) (A) = J\f, K(X\ii)f(n) , (1.4) 

where A is a function of the anisotropy A and the magnetic field h. The kernels are 
given by (as compared to [0 we have performed a similarity transformation, which leaves 
det(l + t/( m )) invariant but changes the kernel V^ m )(A|/i)) 

vim)(xlu) = 1 sinVv) ( el m) (A)el m V) _ e^(X)e^^) ) ( . 

{ ^ } 2tt sinh(A - /i) I (sinh(A - // + 2vq) sinh(// - A + 2irj) J ' 1 ' ' 

= -TV, v (i-e) 

sinh(A — + 2ir]) sinh(A — (jl — 2ir)) 



with functions e±^(A) = (( sinh(A-i^) ) e<? ^ A ' ) ) 2 ■ Here 0(A) is a dual quantum field acting 
in a bosonic Fock space with vacuum |0). The fields are related to canonical Bose fields 
a and at (annihilating the vacuum and its dual respectively, i.e. a(A)|0) = = (0|at(A)) 
via the Bogoliubov transformation 



0(A) = a(A) - / dz/ln 



h(X, v)h(v, A) 



at(iz) , (1.7) 



where h(X,u) = ^^Z^^ 1 ^ ■ The occurrence of the dual fields is a consequence of the 
interacting nature of the Heisenberg model (see Section XI. 1 of [|l|). By construction the 
fields 0(A) commute for all values of spectral parameter [0(A), 0(/^)] = 0. The contribution 
of the dual fields (to the expectation value (0| det(l + V^)\0)) can be visualized by 
decomposing all dual fields according to ( |1.7|) and then moving all exponentials of 
annihilation operators a(l) to the right, picking up contributions whenever passing by an 
at (i/). 

The purpose of the current paper is to derive a system of integro-difference equations 
that drive the Fredholm determinant of the operator t?( m ) in the correlation function 
P(m) ( |1.3|) . This is done in Section |2| and can be considered as an extension of the idea 
of describing quantum correlation function by means of differential equations [Painleve] 
due to E. Barouch, B.M. McCoy, T.T. Wu and C. Tracy and B.M. McCoy @, f§. The 
Fredholm determinant det — ^K^J in the denominator of ( |1.3| ) amounts merely to an 
overall normalization independent of the distance m and will not be considered here. In 
Section ||| the Riemann-Hilbert problem associated with the IDE is formulated. Finally, 
in Section [|, we determine the long-distance asymptotics of the FSFP in some limiting 
cases. 



2 Integro-Difference Equations 

We start by bringing the kernel of V (|1.5| ) to "standard" form Q. We first perform a 
change of variables z = e ^\_~2 (recall that w = e 2l7? ), which maps the real axis on the 
contour C : a — > z = expia (see Fig. |l|) where — ip < a < 2tt + ip (ip < by definition). 
The endpoints £ = e 1 ^ and £ = e~ 1 ^ (we integrate from £ to £) of the contour are related 
to the magnetic field h and the anisotropy A. Using the identity (valid for | < rj < tx 
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and Zi,Z2 G C) 



■ / WZ-I — 1 1 WZn — 1\ 7 

r —Jill) IS b /f, i\ 

as e z i~ w w z 2-™ ; = = = r , (2.1) 

„.,wz\ — l 1 WZ2 — 1 ' V / 



Z\—W W Z2 — W 



the kernel of V^" 1 * 1 is found to be (up to a similarity transform which leaves the determinant 
unchanged) 

(m) / I \ (m) / I \ (m) / i \ (m) , 



V (m) (ziU2) = _± r ds eT\zAs)er\z 2 \s)-e^(z 1 \s)eT\z 2 \s) 
2n Jo z\ — z 2 

where the functions e± l ' ) are given by 

y, , , /.. , . x wz — 1 m ^.A 2 / i W2 — 1 

e!p(zs) = 2" m e"^ ) exp s 

\ z — w ) \w z — w 



e_ 'u s = z(u> — 1 w) z e VK ' ) exp — iw s 

\ z — w J \ z — w 



(2.2) 



(2.3) 



The integral operator t/( m ) now acts on a function f(z) as 

(V™f) ( Zl ) = f dz 2 V^\ Zl \z 2 )f{z 2 ) (2.4) 



where the integration is to be performed along the contour C. We note that V is symmetric 
and nonsingular at Z\ = z 2 . The resolvent R^ of V^ m > is defined by (1+V^ m ')(l— R^) = 
1 and its kernel R^ m \zi\z 2 ) can be written in a form similar to Eq. ( |2.2j) , namely [|TJ 

mz2) = .i. r is i?\'MP\*\')-F\*MfrX*\.) . (2 . 5) 

2tt Jo zi — z 2 

Here f±^ are solutions of the linear integral equations ((1 + V^)/^""^ (z\s) = e^ l \z\s). 



In terms of these functions we introduce the integral operators B^\ a, b — ± acting as 
(BS°/ (m) )(*) = Io°°dt B ( ™\s,t)f(t) with the kernel 

B%\s,t) = ^J c ^ft\z\s)et\z\t), a,b = ±. (2.6) 



The transpose B acts like [B ' fj (s) = J °° dt B(t, s)f(t). We are now in the position to 
formulate the main 
Theorem : 

(i) The lattice logarithmic derivative o/det + V^j is given in terms of the integral 
operator B^ as 

det fl + V (m+1) ) , , u 

\ = det f 1 + B™) (2.7) 
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(ii) The logarithmic derivative of det +y( m )^ respect to the boundaries of the 

contour C is expressed in terms of the functions F^ n \s) = /±"^(£| s ) an d G±(s) = 
/± m) (^|s) as follows 

-zfylndet(l + F) = ^ J" ds - 

+GL(s)fyG + (s) - Gf+(s)^G?_(s)} 

+^ |±| QHds [F + ( S )G_( S )-F_( S )G + ( S )])(2.8) 



(mj TTie following set of completely integrable integro- difference equations for the un- 
knowns F,G,B in (i) and (ii) holds 



(2.9) 



l G (m+D = H G ^-(i-Bt +1) )(B ( r2) T G {m) 



1 ~(m+l) _ 1 jY f . n(m+l) / H (m)\ T \ ^(m) n (m+l) , / D («0\ T \^M 



and 

-i-^B ab (s,t) = ^F a (s)[F b (t) + (F + B^ b -F_B +b )(t)} 

+G a (s) [G b (t) + (G + B_ b - G-B +b ) (*)] J (2.10) 

where a, b — ±. 

TTie additional restrictions necessary to solve the equations uniquely are provided by 
the requirements on analyticity and the asymptotic behaviour of the solutions of the corre- 
sponding linear system (\2.13j) , i.e. by specification of the data in the corresponding 
Riemann- Hilbert problem (see Section^). 

Proof: The proof is analogous to the one for the XXX-case |5| so that we only sketch the 
main steps, (i) is a direct consequence of the shift- equation 

*T* V^\ Zl \z 2 ) z\ = V^( Zl \z 2 ) +-L fda ^\^\s)e M (z 2 \s) 

2n Jo Z\ 



which follows directly from (2.2) and (|2.3|) . (ii) and (iii) follow from the Lax-representation 

1 



1 Am+1) 
~ 


'Ms) 


- f (m) (z\s) - 


1^+1) 


l z \ s ) 


z V + 


and 






—f± 

dip 


(z\s) 


- 1 1 ? 


2tt \z-£ 

+ V 









z 



(m+i) Q x + ( 5 M) T 1 ;M _ ( S M) T /("A (2| S ) , 
1 _ B (m+D] (B^f f^A (z\s) (2.12) 



/±(£|s) / dt (umf-(z\t)-f-mf+(z\t)) 



q.e.d 



It is a remarkable fact that (|2.7|) - (|2.10|) are formally identical to the corresponding ex- 
pressions for the XXX case ||. 



3 The operator-valued Riemann-Hilbert problem 

In this section we show that all results of Section ^| can be reformulated in terms of 
an infinite-dimensional Riemann-Hilbert problem (RHP) for the integral operator valued 
function x( z )- From the solution x( z ) °f tins RHP one can then determine an asymptotic 
expansion of the correlator P{m). 

Let us now consider the following infinite dimensional Riemann-Hilbert problem for 
the integral operator- valued function xi z ) : 

1. x{ z ) is analytic outside the contour C (Fig. [[]). 

2. X~( z ) = X + ( z ) ' L^ mS) (z) for z G C, and x ± are t ne boundary values of the function 
x( z ) as indicated in Fig. [I], and where the (integral-operator valued) "conjugation 
matrix" L^ m \z) is given by 



L (m) (z) = I + 6 m \z) 



(3.1) 



.10. 

where I(s,t) = S(s — t) \ \ and 

1 



i im \z\s,t) 



1 -e { l n \z\s)e { ™\z\t) e { l n \z\s)ef\z\t) 



\ -e {m \z\sy m, (z\t) e ym> {z\s)e K r{z\t) 



(m) 



(m) 



(m) , 



(3.2) 
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3. x( z ) ~ y I as z ~ ¥ 00 • 



In order to simplify our notations we will from now on suppress the m-dependence in 
equations where all quantities are to be taken with the same m and write e.g. £(z\s,t) 
instead of £( m '(z\s,t). In terms of the corresponding kernels the properties 1-3 can be 
rewritten in the following way: 

PI. x( z \ s jt) is an analytic function of z £ C for all s,t. 

POO 

P2. x~(z\s,t) = x + ( z \s,t) + / ds'x + (z\s,s') £(z\s',t) for z G C 

Jo 

P3. x( z \s,t) = I(s,t) + ;*i(s,t) + . . . as z -> oo 

The connection of the RHP 1-3 to the IDE of Section |2] is summarized in the following 
lemma: 

Lemma 1: Suppose now that the solution of the Riemann-Hilbert problem 1-3 exists and 



is unique. Then the junction ^(z) = x( z ) ' 



1 



satisfies the integral operator- 



valued linear system (\2.13i) , (\2.13j) . 



Proof : Applying the standard arguments based on Liouville's theorem and on the m 
independence of the conjugation integral operator L (z) 

I 



we obtain 



Hz) 



Y]/(m+l) f^H 




(3.3) 



-i I \ 1 

• / + - U 
Oil z 



(3.4) 



where U = x (m+1) (0) 



1 




X (0) • As we will now show by determining Uq, 



-i 



(|3.4|) is the integral-operator valued analog of Q2.12|) . 

By construction the operator £(z) (|3.2|) is nilpotent £ 2 (z) = 0, hence we have L~ l (z) 
I — £(z), which in turn implies the equation 



X 1 ( z \s,t) 



Xzi{z\t,s) -Xi2{z\t,s) 
y -X2i(z\t,s) xn(z\t,s) 



(3.5) 
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for the solution of the RH-Problem 1-3. Introducing notations 



x(o) 



B++ 
1 + 



(3.6) 



we conclude from 



that 



x _1 (o) 



l + B T , 
B T _ 



-B 



T 



Bl- 



(3.7) 



Substituting the expansion P3 for x( z ) as z ~^ 00 ( |3 . 4| ) , we see that (Uq)u = 1. 
This together with formulae (|3.6|) and (|3.7|) allow us to rewrite U as 



( 



\ 



B 



(m+l) 



1 



1 



l-B 



(m+l) 



B 



B 



(m) 



^(m+1) Mm) 



T 



(3.1 



We now define functions f±(z\s) via 

( 




dt x + (z\s,t) 




(3.9) 



It can be shown analogously to Section XV. 6 of [|TJ and || that the functions f± defined 
this way are identical to the ones defined in Section |2| for Eq. (|2.5| ). This then implies 
that the kernels B a b(s, t) from ( |3.6|) are just the kernels introduced in Q2.6|) . To prove this 
one has to consider the canonical integral representation for the solution of the problem 

1-3 

dz 2 



1 



2m Jc z 2 — Z\ 
dz 



X + (z 2 ) £(z 2 ) , 



(3.10) 
(3.11) 



which implies 

Taking into account the explicit formula for £(z) and Eq. ( |3.9| ) we obtain the representa- 
tions ( |2.6| ) for B ab (s,t) defined in Q3.6|). This shows that *&(z) fulfills the integral-operator 
valued version of the system (|2.12| ), which can be reobtained from Q3.4|) in the following 
way: rewriting fl3.4|) as 



( -i 

Z 2 



V 



o 



o 

22 




(3.12) 
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and then applying both sides of ( |3.12D to the vector 



V 



z\s) 



and taking into account 



(|3.9|) and that e± 2^ = e± 1+1 ' ) we arrive back at the equations ( |2. 12| ) 



Let us now study the ^-derivative of the function x( z ) m order to connect ( |2.13| ) to 
our RHP. The corresponding analysis is very similar to the ones for the Bose gas |B[ and 
the XXX chain ||. In the neighbourhood of C the function x( z ) can be represented as 



XO) = X{z) ■ Xo{z) 

where x{z) is single-valued, invertible and analytic in that neighbourhood, and 

*>(*)=/- f m^— f £(z) . 



From (|3.13|) we conclude at once that 



dxjz) 
dip 



z-i z-c 



where 



A± = Km (z - z ± ) x(z) ■ ■ Xo\z) ■ r 1 (*) 



Differentiation of (|3.14j ) gives 

( dxo{z] 



Using this together with £ 2 = and ( |3.2|) and (|3~^ ) in (|3.16|) we obtain 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



A ± (s,t) 



z± 
2tt 



(x ■ £ ■r 1 )(z ± \s,t) = ^( X ■ t ■ X- 1 ) (*±\s,t) 



2tt 



/ 



(3.18) 



V 



-f + (z ± \s)f-(z ± \t) f + (z ± \s)U(z ± \t) 
-f-(z±\s)f-(z±\t) f-(z ± \s)f + (z±\t) 
Recalling the definition of the potentials F±(s), G±(s), we can rewrite the final formulae 
for A± as 



A+{s,t) 



J_ -G+(a)G_(t) G + (s)G + (t) 
2tt I -G-(a)G-(t) G-(s)G + (t) 



(3.19) 



4 (s n ^ -JVC*)*'-® i^W^+W 

27T ^ /•• (,)/•• (/) /•• ( S )F+(t) 



From (|3.15|) it follows that *&(z) fulfills the integral-operator valued version of (|2.13|) 

V(z\t',t) . (3.20) 



dip 



Acting with ( |3.20|) on 



urn 



and taking ( |3.19|) into account we arrive at the basic 



equation (|2.13|) for the -^-derivative of f±(z\s). This completes the proof of the Lemma, 
q.e.d. 



4 Asymptotics in some limiting cases 

In this section we obtain the large distance asymptotics of the FSFP correlation functions 
in two limiting cases. The first case corresponds to the limit of very strong magnetic fild. 
In the second case, the magnetic field is arbitrary but the value of parameter 77 is chosen 
to be 37r/4 corresponding to the free fermionic point. 

For very large h close to the critical field h c — 4 cos 2 rj (for which the ground state 
turns into the saturated ferromagnetic state) the integration boundary A in ( |1.4j ) tends 
to zero acoording to 

A = ^\ tan 7]\^Jh c - h + 0(h c - h) . (4.1) 



In the limit h — > ft c , the kernel (|1.5|) can be expanded to order 0{m{h c — ft,) a) as 



V^\m ~ Vb(A|M) = _l ^Hcot,|(A-^)) 

71 X — jl 



(4.2) 



We note that the dual fields do not contribute at all because 

0(O = 0(O) + 0((ft c -ft)i), (4.3) 

and the dual expectation value can simply be dropped. Thus we find that in the limit 
h — ► h c , m — > 00 with m(h c — h) <C 1 Q 

P(m) ~ det (1 + V ) . (4.4) 

The long-distance asymptotics of this Fredholm determinant is known 

1 sin(a(x — y)) 



In det 



1 - 



7T x — y 



~ —-(as) 2 for as — > 00. (4.5) 



1 This inequality rather than m(h c — h)* < 1 is a consequence of a more accurate analysis of the 
contribution of the terms dropped in (|4.2|) to the determinant. 
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In this way we obtain a Gaussian decay of the FSFP for large magnetic fields h « h c at 
large distances 

P(m) ~ e -K^- ft ) m2 for /i -> fr c , (/i c - h)~* < m < (/i c - /i)" 1 . (4.6) 

This result complements the expression for near asymptotics obtained in 

P(m)~e~- v/ram for /i->/i c , m< (h c -h)~^ . (4.7) 

For the case of the XXX chain an analogous analysis can be performed (based on the 
determinant representation ||) giving the same result (|4.6|) (note that h c = 4 for the 
XXX case). 

The above results (and more) can also be obtained from the Riemann-Hilbert problem 
1-3 presented at the beginning of Section |3|. As shown above, in the limit h —>■ h c the 
dual fields can be dropped and the integral kernel ( p..5|) can be approximated as 

| m 

VM(AIai) = - g + (A)g-M -g-(A)e + (/,) } = / riift(A + fr) \ » 

27r sinh(A — /i) \sinh(A — «?y) y 

It is worth emphasizing that this is the only approximation that we will use in this 

approach. Performing the same change of variables as in Sect. ^|this becomes an integral 

operator of the form Q2.4p with kernel 



V^( Zl \z 2 ) = -± ~ e-We + W ^ e±{z) = ^ (4 g) 

2,71 Z\ — Z2 

and the RHP 1-3 above reduces to the ordinary 2x2 matrix RHP without s -integration: 

PI. x( z ) i s an analytic function of z C 

~ ( z~ m \ 

P2. x -(z)= X + {z) 

P3. X (oo) =/. 

The functional determinant of the operator ( [4.9| ) can be rewritten as the determinant of 
anmxm matrix from which the behaviour for small m is found. The large-m asymptotics 
can be obtained in terms of the RHP P1-P3 and yields the following asymptotic formula 
for det(I + V^) |] 



In det (/ + V {m) ) ~ m 2 In 5 for m -> oo 

e* = ^i, (4.10) 
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where A ~ m 2 e . In terms of the magnetic field this condition translates into m{h c — h) 



m 



-2c 



The asymptotic equation ( |4.10|) is actually equivalent to ([Of) as 



In 5 ~ — A 2 cot 2 r/ + 0(A 4 ) , for A -> 

~ -~(h c -h), for h^h c (4.11) 



Equations ( [4.6[ ) and ( [4.1U ) describe the long-distance (m — > 00) asymptotics of P(m) in 



the XXZ model in the limit of small A. 

Let us now consider the free fermionic point 77 = 37r/4 (XXO model) in more detail. 
For this case the RHP P1-P3 is the exact RHP and ( f4.8|) , ( |4.9| ) are the exact integral 
operators for the FSFP for any values of A. Thus fl4.10|) with rj = 3ir/4 is the result for 



the long-distance asymptotics of the functional determinant of / + in the XXO model 
for arbitrary A and thus arbitrary magnetic field h 



In det (i + V {m) ) ~ — In (- + \ — -] for m -> 00 . (4.12) 

V J xxo 2 V2 2cosh2A/ v ; 

For the XXO model we have cosh2A = 2/h. Hence we finally obtain for the asymptotic 
behaviour of the FSFP at the free fermionic point as a function of the magnetic field 

PxxoM ~ f 2 + M 2 (413) 

which coincides with ( |4.6| ) for h —> h c = 2. 

On the basis of the above results we conjecture that the FSFP exhibits a Gaussian 
decay for any value of rj and any value of the magnetic field h < h c . 
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Figure 1 



Figure 1: 

Contour C on the unit circle for the integration with respect to z in the integral operator 
V. '+' and '— ' indicate direction in which the limit z — > C has to be taken to obtain the 
boundary values x ( z ) f° r the Riemann-Hilbert problem in Section |[ 
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